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Abstract 

We compute the shear viscosity at leading order in hot Quantum Electrodynamics. Starting from 
the Kubo relation for shear viscosity, we use diagrammatic methods to write down the appropriate 
integral equations for bosonic and fermionic effective vertices. We also show how Ward identities 
can be used to put constraints on these integral equations. One of our main results is an equation 
relating the kernels of the integral equations with functional derivatives of the full self-energy; it is 
similar to what is obtained with two-particle-irreducible effective action methods. However, since 
we use Ward identities as our starting point, gauge invariance is preserved. Using these constraints 
obtained from Ward identities and also power counting arguments, we select the necessary diagrams 
that must be resummed at leading order. This includes all non-collinear (corresponding to 2 to 2 
scatterings) and collinear (corresponding to 1+N to 2+N collinear scatterings) rungs responsible 
for the Landau-Pomeranchuk-Migdal effect. We also show the equivalence between our integral 
equations obtained from quantum field theory and the linearized Boltzmann equations of Arnold, 
Moore and Yaffe obtained using effective kinetic theory. 
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I. INTRODUCTION 



Shear viscosity is a transport coefficient that characterizes the diffusion of momentum 
transverse to the direction of propagation. It has attracted a lot of attention in the realm 
of relativistic heavy ion collisions lately since it could be an important parameter in Quark- 
Gluon Plasma (QGP) evolution. Its value is still a matter of debate in the heavy ions 

On the theoretical side, computations in strongly coupled 
indicate a low viscosity while kinetic theory calculations 



community (see for example 
Super Yang-Mills theories [2j, 



4| indicate a higher viscosity near the critical temperature in the weak coupling limit. See 
also Ref. 5!] for a comparison of the shear viscosity to entropy ratio in Super Yang-Mills 
and Quantum Chromodynamics (QCD) at weak coupling, showing a significant difference 
between the two near the critical temperature. Recent lattice calculations in pure SU(3) 

1 Q 

6|, 17J show that the shear viscosity is very close to the conjectured Kovtun-Son-Starinet 
bound [3]; however, these studies are not complete since they need some external input to 
solve the ill-posed inverse problem of obtaining the (continuous) spectral density from the 
Euclidean correlator (computed at a discrete set of points). On the phenomeno logical side, 



and transverse momentum correlations 



1 Of ] point 



elliptic flow data (see for example 8 
toward a low viscosity (see also ll| for another interpretation) while others argue that color 
glass condensate initial conditions x3] an d plasma instabilities Jjj] could explain the data 
without invoking a necessarily low viscosity. Viscous hydrodynamic simulations are on the 
way , ll6|, ll7| and will possibly shed some light on these issues. From the above, it is 



clear that understanding the inner workings of shear viscosity (especially in gauge theories) 
is important. 

In the present paper, we are interested in computing the shear viscosity at leading order 
in gauge theories from quantum field theory, starting from the Kubo relation and using 
diagrammatic methods. This calculation has been in done in scalar theory [l8] and has 



been reproduced since then using different methods (e.g. real-time formalism |l9l. l2d. |21| . 
direct ladder summation in Euclidean space 221 ] . 2PI effective action methods 23j and 



quantum kinetic field theory derived from the closed-time-path 2PI effective action 24|). 
The difficulty of the calculation lies in the fact that, due to "pinch" singularities, an infinite 
number of ladder diagrams must be resummed to obtain the result even at leading order. 
The calculation is even more involved in gauge theories; in addition to pinch singularities, 
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the appearance of collinear singularities makes the resummation of another class of ladder 



diagrams necessary 



,4 



25 



261 ] . There exists some attempts at computing the shear viscosity 



in gauge theories from quantum field theory (e.g. real-time formalism 



summation in Euclidean space 



29l | . direct ladder 



2fi|), but to the best of 



22], 2PI effective action methods 
the authors' knowledge, none of these approaches go beyond leading log order accuracy (i.e. 
with corrections suppressed by 0(g ln^ -1 ))) or the large Nf approximation. 

A more convenient way to compute transport coefficients is to use kinetic theory. The 
leading order shear viscosity in hot gauge theories has been computed by Arnold, Moore and 
Yaffe (AMY) using an effective kinetic theory where zero temperature masses are replaced 
by thermal masses {4, 25]. This calculation consistently includes the physics of pinch singu- 
larities, collinear singularities and the Landau- Pomeranchuk-Migdal (LPM) effect. Contrary 
to the scalar case where the equivalence between the quantum field theory approach and 



the kinetic theory approach has been shown 18|, |30j , the landmark calculation of AMY has 



never been verified using quantum field theory. Since the perturbative results of AMY for 
the shear viscosity is often quoted in the RHIC community, it would be good to have a first 
principle quantum field theory proof of this kinetic theory calculation. 

It is thus the goal of this paper to compute the shear viscosity at leading order in hot 
Quantum Electrodynamics (QED) from purely diagrammatic methods and show the equiv- 
alence with the kinetic theory results of AMY. A similar calculation for electric conductivity 
has been performed in 26|, [27j . In the following, we use the same general method outlined in 



Ref. 



26 



271 ]. emphasizing the differences with the electric conductivity calculation. The rest 



of the paper is organized as follows. Section [TT] presents our notation and some background 
material on transport coefficients, both in scalar and gauge theories. We present our con- 
straint on the ladder kernels coming from Ward-like identities in Sect. IIHI Power counting 
arguments are shown in Sect. IIV} in order to determine which rungs should be kept in the 
resummation. The final expressions for shear viscosity, including collinear physics and the 
LPM effect, are presented in Sect. |Vl where we also show the equivalence with the kinetic 
theory results of AMY. We finally conlude in Sect. ED 
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II. BACKGROUND MATERIAL 



A. Notation and definitions 

We present here a summary of our notation along with a list of the various finite temper- 
ature field theory quantities that we use throughout the paper. Latin indices run from 1 to 
3 and represent space components while Greek indices run from to 3 and represent space- 
time components. Boldface, normal and capital letters denote 3-momenta, 4-momenta and 
Euclidean 4-momenta, respectively. We use the metric convention i]^ u = (1,-1,-1,-1). 
The subscripts B, F attached to a quantity refer to its bosonic or fermionic nature (except 
for self-energies and widths, where we use special notations). The subscripts R, I mean real 
or imaginary part and the superscripts ret, adv, cor mean retarded, advanced or autocor- 
relation (i.e. average value of the anti-commutator). A bar over a quantity means that the 
gamma matrix structure is explicitly taken out (e.g. G^(k) = r y fl G{k)). 

Spectral densities are important quantities (especially at finite temperature). Any 2-point 
function can be expressed in terms of them and they satisfy a number of general properties 



[see for example 



31 



32 Q. Spectral densities can be expressed in terms of commutators of 



33 



341]). but it is sufficient for our purposes to give their explicit 



fields (see for example 

expressions in momentum space. Free field spectral densities are given by [351 ]: 

p B (k) = sgn(k°)2n5((k ) 2 - El) (1) 
p F {k) = \27c5{k° - E k )h + (k) + 27r5{k° + E k )hXk)} (2) 

where sgn(/c°) is the sign function, E k = |k|, h±(k) = (7 0= F7- k)/2 and k = k/|k|. Note that 
since we consider systems where the temperature is much larger than any other scale, we put 
m = in the above and all subsequent expressions when the momentum of the excitation is 
hard. From the above expressions (or more generally from CPT invariance), it can be shown 
that the spectral densities satisfy p#(— = — Pb(/c°) and pF^—k) = pF{k) (in the massless 
limit). The delta functions in Eqs. (CD)-© means that the excitations have sharply peaked 
energies and an infinite lifetime. But at finite temperature, any excitation propagating in 
a medium has a finite lifetime due to numerous collisions with on-shell thermal excitations. 
The effect of this finite lifetime is to turn the delta functions in Eqs. (OH)-© into Lorentzians, 
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giving 



l22|: 

pB{k) 
p F (k) 



Ik 



Ik 



•IE, (k° - E k y + ( 7fc /2)2 (hP + E k y + ( 7fc /2)2 

Tfc , /fx , r fc 



;h+(k) + 



(3) 
(4) 



The widths are given by 7fc = ITf^A; = E k )/E k and T k = tr [^^(k = E k )\ /2E k , where 
n(fc) and are the bosonic and fermionic self-energies, respectively. Note that when the 
momentum k is soft, perturbation theory must be re-organized and partial resummation 
of spectral densities is necessary (also called Hard Thermal Loop (HTL) resummations 



36. 



33, 



38]). These resummations give rise to screening thermal masses and may also produce 



Landau damping. In gauge theories, HTLs are also essential to obtain gauge invariant results. 

The analysis in this paper relies heavily on the finite temperature cutting rules. The 
building blocks of these rules are the four propagators of the closed-time-path or "1-2" 
formalism [39, l^. The time-ordered (or "uncut") propagators can be expressed in terms of 



the spectral densities [35j: 
G B /F(k) = 



duo 



(2tt) Pb/f ^ j \k°-uj + ie ^ k° 
where nB/F{k°) are the usual Bose-Einstein or Fermi-Dirac distribution functions. The 



UJ 



1 ± n B/F (u) n B/F (uj) 



uj — it 



(5) 



anti time-ordered propagators are just the comp 



Wightman (or "cut" ) propagators are given by 35] : 



ex conjugate of the time-ordered ones. 



A + B/F (k) = (l±n B/F (k°))p B/F (k) 
A B /F( k ) = ±n B/F (k°)p B/F (k) 



(6) 
(7) 



Our notation for the propagators is different from the one of the "1-2" formalism; the 
correspondence is G = G n , G* = G 22 , A + = G 12 and A~ = G 21 . Using the above 



propagators, we can write down the following cutting rules 



35 



41 



42|]: 



1. Draw all the cut diagrams relevant to the problem considered, where cuts separate the 
unshaded (i.e. "1") and the shaded (i.e. "2") regions. 

2. Use the usual Feynman rules for the unshaded region assigning G B / F (k) to the uncut 
lines. For the shaded region, use the conjugate Feynman rules assigning G* B ^ F (k) to 
the uncut lines. 
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3. If the momentum of a cut line crosses from the unshaded to the shaded region, assign 
A~^, p (k). If the momentum of a cut line crosses from the shaded to the unshaded 
region, assign Ag, F (k). 

4. Divide by the appropriate symmetry factor and multiply by an overall factor of —i. 
These rules are analogous to the zero temperature ones and reduce to them when T = 0. 



Note that the various shadings given by the cutting rules are not a 
related by unitarity (i.e. "vanishing of all circlings" re 



43 



ation 



35 



44 



1 independent: they are 



41 



42) ) and the Kubo- 



451 ] for 3 and 4-point functions 



Martin-Schwinger (KMS) relations (e.g. 3JJ, see also 
KMS relations). 

In some cases, it is convenient to work with linear combinations of shadings. A particu- 
larly useful one is the Keldysh (or r,a) basis in which we can write the physical functions 
(see for example [33| ) : 



■(~ira _ 
tLr B/F = 


iGf /F (k) 


= G B /F{k) 




(8) 


^B/F = 


iG B j F (k) 


= G B/F {k) 




(9) 


IGb/f = 


iG™ /F {h) 




+ &B/ F ( k ) 


(10) 



The G'g'ip is identically zero in the Keldysh basis (due to unitarity). Note that G ra and G ar 
do not depend explicitly on distribution functions. Using the equilibrium expressions for the 
cut propagators (JSJ)- (EJ) , we see that the rr propagator satisfies the fluctuation-dissipation 

[l±2n B/F (k°))p B/F (k) = (l±2n B/F (k )) (iG™ /F (k) - iG% /F (kj) (see 



theorem iG r ^, F ( k) 



for example 



32 



451] ). Note also that any vertex in this basis must involve an odd number of 



a's (see for example 46j), independent of the form of the interaction. To compute a certain 
diagram (topology) in the Keldysh basis, we add all possible r,a configurations consistent 
with the above two constraints. 



B. Shear viscosity from quantum field theory 

Shear viscosity is a transport coefficient that characterizes the diffusion of transverse 
momentum due to collisions in a medium. It is roughly proportional to the mean free path 
of excitations in the medium; simple parametric estimates using kinetic theory shows that 
the shear viscosity in hot QED (i.e. where hot means that the temperature is much larger 
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than the electron mass) behaves as (to leading- log accuracy): 

T 3 



V 



c- 



(11) 



e 4 lne" 1 

where T is the temperature of the medium, e is the electromagnetic coupling constant and 
C is a numerical coefficient that can only be obtained from a detailed analysis {4, [25, \^\ . 
The coupling constant dependence of Eq. (TlTI) is expected, since transport coefficients are 
roughly proportional to the mean free path and thus inversely proportional to the scattering 
cross section of the processes responsible for transport (Coulomb scattering in the present 
case) . 

From linear response theory and Ward identities, it is possible to express transport coef- 
ficients in terms long distance correlators between conserved currents. The resulting Kubo 



formula for the shear viscosity rj is given by 



IS 



48 



49|: 



77 = — lim / <fxe ik - x (vr,- 7 (t,x x 
' 20 k=o ' 



7T 



'J 



(0)) 



eq 



(12) 



where ^(x) 



Tij(x) — 5ijT-(x)/3 is the traceless part of the stress tensor. Note that the 



averages in Eq. (fl2l) are done with respect to an equilibrium density matrix even if the shear 
viscosity is a non-equilibrium quantity; Eq. ffl2|) thus lends itself to a diagrammatic analysis 
and allows one to compute the shear viscosity from first principles. 



In QED, the (traceless) stress tensor has the form (at leading order in the coupling) 22]: 



TT ij {k) 



25 ij 



( 7 -k) 



k l k? - 



k 2 5 ij 



5 st A t 



(13) 



where ip and A s are electron and photon field operators, respectively. From Eq. ( fl3l) . we see 
that such an operator insertion leads to both bosonic and fermionic vertices. Figure [U shows 
the diagrammatic expansion of the Kubo relation ( Tl2l) . Naively, only one- loop diagrams 
should be kept at lowest order, implying no mixing between bosonic and fermionic vertices. 
But these diagrams suffer from so-called "pinch" singularities that make their evaluation 
more complicated than it appears. The argument goes as follows. The low frequency limit in 
the Kubo relation (fl2l) gives rise to products of propagators G(p) with the same momentum. 
Due to the pole structure of finite temperature propagators (one in each quadrant), one faces 
situations when the integration contour is "pinched" between two poles on opposite sides of 
the real axis in the complex p° plane when two propagators with the same momentum are 
multiplied together. In such a case, viscosity diverges as i] ~ J dp G{p)G{p) ~ 1/e, where e 



p+q 




P 

+ — f rv- + (-) 



p+q 



+ U v q + V + (•••) 



p 

+ -CD- + " 

FIG. 1: Diagrammatic expansion in ladder diagrams of the Kubo relation (|12j) . Dotted lines 
represent stress insertions and grey squares represent 4-point functions (called "rungs" ) . When the 
external momentum q goes to zero, the two "side rail" propagators have the same momentum and 
produce a "pinch" singularity. 

is the usual infinitesimal ie prescription. These "pinch" singularities are of course artificial: 
since at finite temperature excitations always suffer collisions with thermal particles coming 
from the medium, they have a finite lifetime (i.e. ie is effectively replaced by iT). This 
finite lifetime regulates the pinch singularities and makes the shear viscosity proportional to 
1/r ~ l/e 4 me _1 , similar to the parametric estimate (TTTT) . 

The use of resummed propagators to regulate pinch singularities also makes the resum- 
mation of an infinite number of ladder diagrams necessary to obtain the correct leading 
order result. One way to do this infinite resummation is to write the Kubo formula in terms 
of two effective vertices (bosonic and fermionic) that satisfy two coupled integral equations. 
Schematically, we have: 



V 20 



/ ^ I F(k)^(k)V F (k) + j ^-T B {k)T B {k)V B {k) 
with the accompanying integral equations: 

V F (k) = I F {k) + j JC (1) (k,p)V F (p) +/ £<2)(M2>*(P) (15) 



d A p . . , . f d A p 



V B {k) = l B (k) + Jj^z IC {3) (k,p)V F (p) + J j^i JC {i) (k,p)V B (p) (16) 

where /C@) = M.(i)J-. See Fig.[2]for a graphical representation of Eqs. (fT4|) - (|T6|) . The symbols 
Fb/f represent pairs of bosonic/fermionic "side rail" propagators (note that the "ladder" 
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where 





If M (1) F f D f M (2) F b D b 



+ 






I B M (3) F F D F 



+ 




M (4) F B D B 



^i — r 




FIG. 2: Schematic representation of Eqs. ()14|) -(fT6|). The symbols refer directly to the equations: 
Ib/f are stress tensor insertions, M.^ are rungs, Tbif represent pairs of side rail propagators and 
T>B/p are effective vertices. 

diagrams in Fig. [2] are on their sides, meaning that the side rails are on the top and bottom 
of the diagram). In the limit q — > 0, the two side rail propagators that connect M. and D 
have the same momentum and pinch, producing a 1 /T p or a l/j p factor. The "rungs" M. are 
amputated 4-point functions that must be of the same order as T~ x but otherwise arbitrary. 
Note that there are four categories of rungs that allows the mixing of bosonic and fermionic 
effective vertices. The bare vertices Xb/f represent amputated stress tensor insertions. The 
effective vertices T>b/f encode the information about the infinite resummation of ladder 
diagrams. Closing them with the appropriate bare vertex and side rails, we obtain the Kubo 
relation (TUEj) . 

The nonperturbative resummation of an infinite number of ladder diagrams due to pinch 
singularities is a common feature of both scalar and gauge theories. However, as is usu- 
ally the case, the gauge theory calculation contains other complications, among them HTL 



resummations and gauge invariance. These issues are discussed in Ref. [26| in the case 
of electric conductivity in QED; we briefly come back to the issue of gauge invariance in 

Sect, nn 

Another very important complication of gauge theories is the sensitivity of transport coef- 
ficients to soft and collinear physics even at leading order 4|, l25|]. As with pinch singularities, 



collinear singularities generically appear when two propagators with momenta p and p + q 



9 



that are nearly collinear (i.e. p • q ~ 0(e 2 T 2 )) are multiplied together and integrated over: 



dp" G(p)G(p + q) ~ / dp" G IC \p)G adv (p + q) 

dp ' ' 



(p° + ie) - E P J V (p° + g° - ie) - E 1 



p+q 



(17) 



g° + (E p+q - E p ) - 2ie 

There are two ways to get a divergence in Eq. (ITT)) . First, in the limit q — > 0, the ex- 
pression diverges as 1/e and corresponds to the pinch singularity case discussed previously. 
Second, when q is nonzero but nearly on-shell and the angle between the quasiparticles is 
parametrically small (i.e. 9 pq ~ 0(e)), we have E p+q ph E p ± |q| and the expression diverges 
as l/0 pq (or 1/e in the perfectly collinear case). The regularization of this new singularity 
using resummed propagators introduces coupling constants in the denominator and changes 
the power counting dramatically. In the scalar case where there are only pinch singularities, 
only one resummation of an infinite number of ladder diagrams with a finite number of 
rung types is necessary [181 ] . In gauge theories where two types of singularities are present 
(pinch and collinear), two resummations are necessary. The one due to pinch singularities 
is similar to the scalar case. The second resummation takes place inside particular types of 
rung where collinear singularities are present, making the number of rung types effectively 
infinite. This infinite class of rung types is the manifestation of the LPM effect in our ap- 
proach. In practice, the two resummations are done using integral equations embedded in 
each other: the integral equations for pinch singularities are shown in Eqs. (115j) - (116j) and 
the ones for collinear singularities is hidden in the four kernels K.^ . We come back to these 
issues and the subtleties of power counting in gauge theories in Sect. [IVj 



III. CONSTRAINTS ON THE LADDER KERNELS FROM WARD-LIKE IDEN- 
TITIES 

Ward identities relate N-point correlation functions to a linear combination of lower- 
point functions and are a direct result of gauge symmetry. In QED, Ward identities must be 
respected in order to preserve transversality (and thus gauge invariance). A particular case 
of interest to us is the Ward identity relating a vertex to a combination of two propagators. 
Since the resummation of ladder diagrams in transport coefficient calculations are written 
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FIG. 3: Momentum convention for the Ward-like identities (| 18 |) - f)19|) . Dotted lines stand for energy- 
momentum tensor insertions and blobs represent the amputated effective vertices T> (bosonic or 
fermionic) . 



in terms of effective vertices satisfying integral equations (c.f. Eqs. (TToT) - (TToT) ). we expect 



Ward identities to put constraints on these integral equations. It is shown in Ref. [26| how 
to obtain these constraints for electric conductivity in QED. In this section, we sketch how 
to obtain such constraints for shear viscosity in QED; we refer the reader to Ref. [26] for 
details (no conceptual difficulties are added in the shear viscosity case). 

In general, Kubo relations express transport coefficients in terms of different conserved 
current insertions. For electric conductivity, the conserved current is the usual electric 
current j M a nd the corresponding Ward identity (relevant for transport calculations) can be 



derived 



26 



501 ]. For shear viscosity, the conserved current is the energy- momentum tensor. 
Using standard procedures (e.g. [5l|), we can similarly derive two Ward-like identities for 
a T^ v insertion, one for each type of effective vertex. In Euclidean space, the result is (see 
Fig. [3] for the momentum convention): 

Q U V F U (K,K-Q) = (K»-Q»)G F \K)-K»G F \K-Q) (18) 
Q»V%ap(K, K-Q) = {K»- Q»)G B l aP {K) - K»G B \ p {K - Q) (19) 

where T>b/f are the bosonic/fermionic amputated effective vertices defined in Sect. Ill Bl 
To go from Euclidean space to Minkowski space, we need to analytically continue K and 
Q towards real energies. The proper choice here is dictated by the physics of transport 
coefficients encoded in the Kubo relation. First, the Wightman correlator in Eq. (fT2l) can 
be expressed as the imaginary part of a retarded correlator. Second, in the pinch limit, 
the two side rail propagators that connect the ladder kernel to the effective vertex have 
different boundary conditions, namely G r F t {k)G F ^ v {k — q) or G F hr {k)G r F t {k — q). These two 
requirements uniquely fix the analytic continuation to K — > k° + ie, K — Q — > k° — q° — ie and 
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Q — > q° + 2ze. Taking the pinch limit (g° — > 0, q — > 0) after doing the analytic continuation 
and using G F \ ct/adv (p) = [ 7 °(p° ±*T p /2) -7 ■ p] and GT^^ip) = ^[(pP ±i 7p /2) 2 - |p| 2 ] 
(valid when p is nearly on-shell), Eqs. (ITS1) - (ITU1) become 

lim q v V F v (k + ie,k-q-ie) = ik^°T k = 2ik^f (k) (20) 
lim q v V%Jk + ie,h-q-ie) = 2%k»g af3lk = 2iFIT;^(k) (21) 

q^O ' 

The last equality is valid near k° pa |k|. These last equations relate the effective vertex of 
the integral equation to the imaginary part of the on-shell retarded self-energy in the limit 
relevant to transport in the case of an energy-momentum tensor insertion. Note that we 
used fully resummed retarted/advanced propagators because of the need to regularize pinch 
singularities. 

Let us now see what are the implications of the Ward-like identities (l2Up - (l21l) on the 
integral equations f|T^]) - ffTB]) . Starting from the Euclidean versions of Eqs. ([15]) - (}!?)]) . 

V^(K,K-Q) = ZJ»(K,K-Q) + J (K, P, Q)V?{P, P - Q) 

+ Sw) i *$ ( *' P ' Q)v ^ {p > p - Q) (22) 

V% af3 (K,K-Q) = T^ a/3 (K,K-Q) + J ^JC maP (K,P,Q)V^(P,P-Q) 

+ J^ K w( K > P >® V Z<#( P > p -® (23) 

we multiply both sides by Q u and do the analytic continuation that leads to pinch singulari- 
ties. Note that this step is delicate, since the sum over Matsubara frequencies must be done 
before the analytic continuation; the result is that the integral equation keeps its form in 
Minkowski space, as shown in [26( for electrical conductivity. After taking the q — > limit, 
we use the Ward-like identities (I2~D]) -([2T1) on both sides of Eqs. (122>(l23"]) . We thus obtain: 

2*"£f(k) = J ^ M (1) (k,p) p^pUp) +/(^4 M<*(k,p) P»P + BaP ( P ) (24) 

2Pn^(k) = |^^ w ^(M^(p)+/^^(4)(M^p^W (25) 

where we have used the decomposition fC(k,p,q) = M.(k,p,q)G(p)G(p — q) and Pg/ F (p) = 
i(G V g, F {p) — G^J F {p)). Here S} et (k) and IPpJ^k) are the imaginary parts of full retarded 
self-energies and the .M^'s represent coherent additions of rungs. These last equations 
implement the constraints imposed by the Ward-like identities (|20p - (j2~T]) ; since they are 
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expressed only in terms of known quantities (self-energies and spectral densities), they also 



impose constraints on the Mia's. Following the procedure in [26j, it is possible to invert 
the system of equations (1241) - (1251) and express the Ai^s in terms of functional derivatives 
of self-energies. The result is (Lorentz indices for bosonic self-energies and spectral densities 
are not shown for simplicity): 

r{1 + b)MmHM = ^{^^-) ^^fc^p^M) (26) 

^+^<^>^Nir) -^'^Nlf 1 ) (27) 
''" +i ) w «*')=''(w ! ) ^^'ft^K^w 1 ) (28) 

+ f)M m „,(k, q) = fc» ( ' 5(2 5 "^ ( ) <;)) ) -» ^)(*.9)«* y ( j ^P ) ( 29 ) 

where the rung kernels and self-energies are expanded in terms of their number of bosonic 
(b) or fermionic (/) spectral densities: 

oo oo oo 

m = M bf sf = vfif nf = J] nyg, (30) 

b,f=l bj=l bj=l 

These definitions mean that Ai, Sj et and U} et are "blobs" containing diagrams of all orders; 
these diagrams can be reorganized in terms of their number of bosonic or fermionic spectral 
densities. The proportionality relations are statements about the structure of the Ai^s. 

The diagrammatic interpretation of Eqs. fl26|) - fr29|) is quite natural. Since the propa- 
gators (j31)-([7)) are all proportional to spectral densities, then functional derivatives with 
respect to Pb/f can be interpreted as opening boson/fermion lines in a Feynman diagram. 
Thus each rung kernel is obtained by opening bosonic/fermionic lines in the appropriate 
self-energy diagrams. Equations fT2oT) - fl2^1) are the desired constraints on the rung kernels. 

The above constraints implement the physics of transport coefficients, namely the ap- 
pearance of pinch singularities due to the low frequency, low momentum limit. They serve 
two main useful purposes. First, they serve as a guide to find the necessary rungs in order 
to obtain the shear viscosity at the desired level of accuracy. The recipe is simple: make a 
loop expansion of all the bosonic/fermionic self-energies of the theory, keep the self-energies 
up to the desired level of accuracy and then open them in all possible ways to get the 
rungs. This recipe seems to work well in theories where only pinch singularities are present 
(such as scalar theories), although it is hard to say if it works at all orders, since there 
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exists no transport coefficient calculations that go beyond leading order (except for large N 



theories 



52j]). When other singularities are present (such as collinear singularities), some 



rung resummation is necessary and a loop expansion is thus not sufficient; additional power 
counting arguments must be supplied in order to get leading order results. We come back 
to these issues in more details in Sect. HVl 

A constraints similar to Eqs. (|26|) - (|29|) has been obtained in the case of electric conduc- 



tivity in QED 



26]. In that context, the starting point to obtain the constraint is the usual 



Ward identity for a current insertion. Since the Ward identity for a current insertion is 
a direct consequence of gauge symmetry, the resulting constraint is an implementation of 
gauge invariance for the electric conductivity calculation. More precisely, the constraint tells 
us which self-energies must be resummed in the side rail propagators for each rung present 
in the integral equation kernel so as to keep everyting transverse. In the shear viscosity case, 
the constraints (126|) - (|29|) are obtained from spacetime symmetries, not gauge symmetries. 
Thus the interpretation of these constraints as an implementation of gauge invariance is not 
direct. We argue that the constraints ( I26l) - (l29l) are sufficient to preserve gauge invariance. 
A quick argument would be that the constraints ( I26l) -( l29l) have the same form as the one 



for electrical conductivity 26J, so they should play the same role. We show in more details 
in Appendix [A] that the shear viscosity is indeed gauge parameter independent. 



We note t 



rat expressions similar to Eqs. (T26|) - (T29|) are obtained using 2PI effective actions 



methods [23, 



281 ] , with notable differences. In 2PI methods, the "constraint" is in coordinate 



space and comes naturally from standard functional relations. The important point is that 
the kernel of the integral equation is given by the functional derivative of the self-energy with 
respect to a dressed 2-point function, where the self-energy is itself given by the functional 
derivative of all amputated 2PI diagrams with respect to dressed 2-point functions. No 
reference is made to the low frequency, low momentum limit or pinch singularities. We also 



mention that there seems to be gauge invariance issues with 2PI methods because o 



' the need 



to truncate the 2PI effective action 



53 



54 



551 ] (see also the developments in [56 



5jj). In 



contrast, our method does not make any reference to 2PI effective actions and start directly 
from symmetry principles, the Ward identity being the expression of these symmetries for 
quantum mechanical amplitudes. It is more specific in the sense that the physics of transport 
coefficients is necessary here to obtain our constraint. 
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IV. POWER COUNTING 



In this section, we present power counting arguments to supplement the constraints (126]) - 
|2"9|) . We separate the analysis in two parts. In the first part, we follow the recipe used in 
26| and outlined in Sect. IIH] i.e. we expand the self-energies of the theory, keep only the 
leading order ones and open them according to Eqs. (j26!) - (!29l) . We then use power counting 
arguments to verify that the rungs obtained contribute at leading order. In the second 
part, we show that the naive expansion in terms of coupling constants is not sufficient when 
collinear singularities are present. As pointed out in Sect. Ill B\ a restricted but infinite class 
of diagrams must be resummed in order to get leading order results. We again use power 
counting arguments to identify this class of diagrams. For the rest of this paper, we divide 
the rungs in two categories, the ones containing collinear singularities (J\f) and those that 
do not (M). 



A. Power counting without collinear singularity (A4) 

We first consider the case where the rungs do not contain collinear singularities. As can be 
seen from Eqs. ( !T5|) - (|T6|) . there are four different types of rung kernels Aiu). We treat them 
separately in the following. According to Eq. ( [26]) . the rung kernel A4n\ is given by opening 
electron propagators in the imaginary part of the retarded self-energy of the electron. At 
one-loop, the imaginary part of the electron self-energy is zero since an on-shell massless 
excitation cannot decay into two on-shell massless excitations; it is thus necessary to go to 
two loops for a leading order result (note that one-loop imaginary self-energies are possible 
when the photon is exactly collinear to the electron, a case we consider in Sect. HVBj) . The 
two-loop imaginary self-energies of the electron with their corresponding rungs are shown 
in Fig. HI Note that there are many more cuts that correspond to imaginary two-loop self- 
energies. We consider these other cuts when writing down the integral equation for the 
effective vertices (c.f. Eqs. ( |36|) - (l37|) ); for the moment, we are only interested in the rung 
topology. The power counting for the rungs in Fig. H]is done in Ref. |26(] . The result is that 
the rungs corresponding to Coulomb scattering ((al),(bl),(cl)) are 0(e 2 ) and the others 
((dl),(gl),(jl),(kl),(ll)) are 0(e 4 ) when the off-shell exchange momentum is soft; but due 
to a partial cancellation between 2 — > 2 processes that do not change species types, all rungs 
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FIG. 4: Rungs of Aim, obtained by functionally differentiating with respect to a fermionic spectral 
density the electron self-energies in the left column (c.f. Eq. (|26p ). Power counting arguments (see 
Ref. [26f] ) show that all these rungs contribute at leading order except diagrams (el) and (fl), 
because massless three-body on-shell decays are suppressed. Note that rungs (hl)-(il) and (kl)- 
(11) are different cuts of the same topology; any cut can be used for power counting purposes 
because, in the end, all possible cuts of a given topology must be included in the integral equation 
(c.f. Eqs. (f36j) - ([3T|) ). Since cut propagators represent on-shell excitations, all rungs (except rungs 
(el),(fl),(hl),(il)) are equivalent to 2 — > 2 scatterings with a soft exchange. Note also that in 
a certain kinematical regime, rung topologies (gl)-(il) contain collinear singularities and can be 
converted into 1 — > 2 collinear scatterings (we come back to this issue in Sects. IIVBI and IV Bp . 



Figure taken from 



2m- 



are effectively 0(e 4 ). 

The rung kernel AAm is obtained by opening photon propagators in the same electron 
self-energies used for Aim- The resulting rungs are shown in Fig. [5j Let's estimate the size 
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of rung (c2), reproduced with momentum labels in Fig. [6j The expression for the rung is: 

M {c2) = -*e 4 f ^A B (k - l)A B (l - p)G* F (l)G F (l) (31) 

where for clarity we omitted Dirac matrices and Lorentz indices (irrelevant for noncollinear 
power counting). If all momenta are hard, then (c2) is 0(e 4 ) (this is true of all the rungs 
considered here). The size of the rung is the same when the loop momentum is soft. To 
see this, consider the momentum regime where k, p are hard and on-shell while I is soft 
and off-shell. Using the delta functions inside the cut propagators to do two integrals over 
I, we get a d 2 l ~ 0(e 2 T 2 ) suppression from phase space. This phase space suppression is 
compensated by the two fermionic propagators Gp(l), which are both 0(e~ 1 T~ l ) in size 
when I is soft. Rung (c2) is thus 0(e 4 ) for both hard and soft loop momentum. The size of 
rung (d2) is obtained similarly and is also 0(e 4 ). 

The power counting of rungs (e2) and (f2) is done in the same way. For example, let's 
consider rung (e2). The expression for the rung is (the momentum labels are shown in 
Fig.!]): 

/ -§tf M w = ' ie ' SwfS (l^ A * (/)A ^ + 1 ~ l)G * Ak " p)Gp{k + 1) (32) 

where we explicitly write the integration over p coming from the integral equation (c.f. 
Eqs. ( fl5|) - (|T6|) ). In the momentum range where k, p, I are hard and on-shell (this last condi- 
tion being automatically enforced by the delta function in A B (l)), the fermionic propagators 
are generically 0(T~ 2 ). With the additional requirements that k ■ p ~ k ■ I ~ 0(e 2 T 2 ), we 
have that (k — p) 2 ~ (k + I) 2 ~ 0(e 2 T 2 ) and both propagators are 0(e~ 2 T~ 2 ). On the 
other hand, these additional requirements are equivalent to 6k P ~ &ki ~ 0(e) (where 8 a b 
is the angle between momenta a and b), implying that phase space is restricted such that 
d 3 p ~ |p| 2 sin 8kpd\p\d8k p d<f) ~ 0(e 2 T 3 ) (and similarly for d 3 l). Combining all factors, we 
see that rung (e2) is 0(e 4 ) in all momentum regimes where the exchange propagators are 
off-shell. Note that this topology also contains collinear singularities in a certain kinematical 
regime; the consequences of this fact are explored in Sect. IIVBI 

The analysis of the rung kernel Ai(s) is similar to Ai(2)- Figure [7] shows the rungs 
contained in Ai^) obtained by opening fermion propagators in the imaginary part of the re- 
tarded self-energy of the photon (see Eq. (|28|) ). To do the power counting of -M( 3 ), notice that 
the rungs in 7W(2) and .M(3) are "mirror images" of each other, i.e. the incoming/outgoing 
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FIG. 5: Rungs of Mr 2)1 obtained by functionally differentiating with respect to a bosonic spectral 
density the electron self-energies in the left column (c.f. Eq. (j27[) ) . Note that (e2)-(g2) and (f2)-(h2) 
are different cuts of the same topologies; any cut can be used for power counting purposes because, 
in the end, all possible cuts of a given topology must be included in the integral equation (c.f. 
Eqs. (I36p -(|37p), Except for rungs (a2) and (b2) containing suppressed three-body decays, they all 
contribute at leading order for shear viscosity and correspond to 2 — > 2 scattering processes. Note 
that in a certain kinematical regime, rung topologies (e2)-(h2) contain collinear singularities and 
can be converted into 1 — > 2 collinear scattering processes (we come back to this issue in Sects. IIVBI 
and|VBD. 

electrons on the left of the rungs in Mm become outgoing/incoming electrons on the right 
in .M(3) (similarly for the photons). Exception to this rule are rungs (a2)-(b2), (b3)-(c3), 
(f3)-(g3); they are however suppressed by three-body decays and are not included in the 
analysis. Note also that A4(g\ contains twice as many rungs compared to Mm- This is due 
to the fact that both electrons and positrons can flow in the RHS fermion loop in A4(3), 
whereas it is only a photon loop in M.m. Coming back to the power counting, we can see 
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FIG. 6: Momentum labels used to do the power counting of rungs (c2) and (e2) (c.f. Fig. [5]). 
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FIG. 7: Rungs of -M( 3 ), obtained by functionally differentiating with respect to a fermionic spectral 
density the bosonic self-energies in the left column (c.f. Eq. (|28p ). Except for rungs (b3)-(c3) and 
(f3)-(g3) containing suppressed three-body decays, they all contribute at leading order for shear 
viscosity and correspond to 2 — > 2 scattering processes. Note that in a certain kinematical regime, 
rung topologies (i3)-(p3) contain collinear singularities and can be converted into 1 — > 2 collinear 
scattering processes (we come back to this issue in Sects. ITVBl and IV Bp . 



that all the power counting done for Ai(2) can be applied to A4(3) since they both share the 
same "mirror image" rungs. 

Figure [8] shows the rungs in A4u), obtained by opening photon propagators in the imag- 
inary part of the retarded self-energy of the photon (see Eq. (129~1) ). The expressions for 
the rungs (a4)-(d4) are almost identical to the ones of (c2) and (e2); the power counting 
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FIG. 8: Rungs of Mu\, obtained by functionally differentiating with respect to a bosonic spectral 
density the bosonic self-energies in the left column (c.f. Eq. (|29p ). All the rungs contribute at 
leading order for shear viscosity and correspond to 2 — * 2 scattering processes. Note that none of 
these rungs enter the (leading order) collinear analysis of Sect. IIVBI 

is thus very similar and gives 0(e 4 ) for all four rungs. This is to be expected, since there 
are no Coulomb divergence (i.e. no soft boson exchange) in (a4)-(d4) that would make 
the rungs 0(e 2 ). Note that none of the rungs in Fig. [8] are included in the (leading order) 
collinear analysis of Sect. IIVBI since they all pick up extra factors of coupling constants in 
the collinear regime. We come back on this issue in the next section. 



B. Power counting with collinear singularities (J\f) 

The power counting with collinear singularities is more subtle than the one done in the 
previous section, due to the additional restrictions to a particular kinematical regime. The 
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following discussion is largely based on the work of Arnold, Moore and Yaffe (in particular 



Ref. 



58|) 



In Ref. |26(, we use the "1-2" formalism to estimate the size of all rungs. However, the 
1-2 formalism is not the most appropriate tool for power counting, especially when collinear 
singularities are involved. For instance, all four "1-2" propagators have a temperature de- 
pendent part that can dominate in certain kinematical regimes. In comparison, the Keldysh 
basis (c.f. Eq. [8]) is cleaner, since only the G rr contains distribution functions; it is thus 
easier to spot Bose-Einstein enhancements when momenta are soft and evaluate the size of 
each propagator. The Keldysh formalism has another advantage compared to the "1-2" for- 
malism. The fact that the Keldysh basis is simpler (i.e. absence of aa propagators, absence 
of vertices with an odd number of a's) combined with the easiness of identifying pinching 
contributions (i.e. G ra G ar or G ar G ra ) makes it very convenient. In particular, it renders 
apparent some near cancellations that are otherwise very hard to see in the "1-2" basis. We 
thus use the Keldysh basis to identify leading order rung topologies in the following; once 
one is found (and there are no near cancellations), we can include it (and all its cuts) in the 
final integral equation. 

The first step is to identify the leading order fermionic and bosonic self-energies when the 
initial hard electron or photon is in the collinear regime. As explained in Sect. Ill B\ collinear 
singularities occur when a retarded and an advanced propagator with nearly collinear mo- 
menta are multiplied together. For power counting purposes, we also need at least one of 
the photon propagator to be rr (we come back to this point below). By inspection, we see 
that only the self-energies shown in Figs. l9"lfT0l satisfy these criteria (note that the photon 



self-energies in Fig [10] are identical to the ones of 58(). The rest of this section is devoted 
to the power counting of the resulting rungs. 

The first class of rungs comes from opening fermionic self-energies (see Fig. [9]). Rung 
(ul) is reproduced with r,a indices and momentum labels in Fig. [TT] (rung (xl) is done in a 
similar way). The expression corresponding to rung (ul) is: 

xi a G r F a (p + l)^V F {p)^G^{k - l)Y (33) 

where we explicitly write the integral over p and the effective vertex coming from the integral 
equation (c.f. Eq. [15]). Note that there are other possible r,a assigments (this is why we put 
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FIG. 9: Rungs of A/i and A/2 obtained from the fermionic self-energies in the left column by 
functional differentiation with respect to pp and ps, respectively. A thick line corresponds to a 
soft propagator. From power counting arguments, it can be shown that only rungs (ul), (xl), (u2), 
(v2) contribute at leading order and must be included in the collinear analysis. 



N (31 (u3) (v3) (w3) (x3) 




FIG. 10: Rungs of A/3 and A/4 obtained from the bosonic self-energies in the left column by 
functional differentiation with respect to pp and ps, respectively. A thick line corresponds to 
a soft propagator. From power counting arguments, it can be shown that only rungs (u3)-(x3) 
contribute at leading order and must be included in the collinear analysis. 

a proportionality sign in Eq. fl33l) ): some of them are sub leading, some of them contribute 
at leading order, but the important point is that no near cancellations between different 
r,a configurations are at work for this rung. The external legs are hard and nearly on-shell 
due to pinch singularities, i.e. k ~ p ~ T and k 2 ~ p 2 ~ 0(e 2 T 2 ). We consider the 
case where the loop momentum I is soft (the case where I is hard corresponds to the non- 
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FIG. 11: Momentum and r,a labels used to do the power counting of rungs (ul) (c.f. Fig. [9] and 
Eq. (|33|) ). A thick line corresponds to a soft propagator. 

collinear rung (gl)). The integral over dl° is dominated by the kinematical range where both 
fermionic propagators are collinearly singular; this happens over the parametrically small 



frequency width dl° ~ 0(e 2 T) [58| and when (p + l) 2 ~ (k — l) 2 ~ 0(e 2 T 2 ) (or, equivalently, 
dpi ~ dki ~ 0(e)). These last requirements are equivalent to the statement that the side rail 
electrons must be collinear with the hard photon. In particular, Q v \ ~ 0(e) implies that the 
phase space of p is restricted to an 0(e 2 ) region since d 3 p ~ |p| 2 sin p /(i|p| d9 p id<fi ~ 0(e 2 T 3 ). 

When / is soft, G r ^(l) = (1 + 2ub(1°))pb(1) is HTL resummed and its size depends on 
the momentum flowing through it. The dominant contribution comes from soft spacelike 
momenta. In this situation, Landau damping gives rise to an 0(e 2 ) imaginary self-energy; 
the size of the spectral density is thus Pbhtl(0 ~ nr(Z)/(Z 2 + IL-(Z)) 2 ~ 0(e~ 2 ) and the 
size of the propagator is G'b HTL (Z) ~ 0(e~ 3 ). Collecting all powers of e, we get e 4 x e 2 x 
(e 2 x e 3 ) x e~ 3 x (e~ 2 x e -2 ), where the e 4 comes from the four explicit vertices, e 2 from 
the (small angle) restriction on the phase space of p, (e 2 x e 3 ) from the (parametrically 
small) integration over 1° and the soft integration over 1, e~ 3 from the soft HTL resummed 
rr propagator and (e~ 2 x e~ 2 ) from the two collinearly singular fermionic propagators. Rung 
(ul) is thus 0(e A ) and contribute at leading order in this particular regime. 

Rungs (yl)-(zl) and (vl)-(wl) have the same topology but are evaluated in different 
momentum regimes. We argue that the momentum regime of rungs (yl)-(wl) (soft momen- 
tum exchange between the side rails) is suppressed compared to the case when the side rail 
momentum exchange is hard. A quick argument to see this suppression is to notice that for 
any rung with incoming momentum k and a soft exchange (p ~ 0(eT)) between the side 
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rails, we can expand the two pinching side rail propagators that connect the rung to the 
effective vertex in a Taylor series, i.e. G(k + p) ~ G(k) + p' J 'd II G(k). The first term in the 
Taylor expansion gives rise to a rung with incoming and outgoing momenta k and so does 
not disturb the ladder at all. We are thus left with only the term proportional to p that is 
0(e) suppressed. In other words, the change induced by adding a rung with a soft exchange 
between the side rails is parametrically smaller than the one caused by adding a rung with a 
hard exchange. This argument is valid for all rungs with no hard exchange between the two 
side rails. A better way to see this suppression is to look directly at the resulting linearized 
Boltzmann equation and see the cancellation between scattering processes for which both 



incident excitations undergo a soft scattering without changing species type [26l. |47|. 



The complete power counting of rungs (vl) and (wl) is done in Ref. 26|]. The result is 
they are both 0(e A ). They should nevertheless not be included in a leading order analysis 
due to a near cancellation 27j 63j. The cancellation can be seen as follows. The assigment of 
r,a labels in rung (vl) is tightly constrained because of the need for the side rail propagators 
(connecting the vertex to the rung) to pinch, the need for the two fermion propagators inside 
the vertex correction to pinch and also the need for the soft bosonic vertex correction to be 
rr. The two possible r,a configurations of the "rung-side rails- vertex" system are shown in 
Fig. [12j Adding the two contributions, we get: 

M {vl) ~ e 4 (Common part) [G r F r (p + l)G a ;(p) + G r F a (p + l)G r ;(p)] 
~ e 4 (Common part) [(1 - 2n F (p° + l°))G r ^(p + l)G a F r {p) 

-(l-2n F (p°))G r F a (p + l)G a ;(p)] (34) 

where, to obtain the second line, we used G r F (p) = (1 — 2n F (p )) (G r F (p) — G c f(p)) and kept 
only the pinching part. When I is soft, the two contributions nearly cancel (up to 0(e) 
corrections) and make the rung subleading. 

Rung (u2) is reproduced with momentum and r,a labels in Fig. [121 The corresponding 
expression is: 

(35) 

where we explicitly write the integral over p coming from the integral equation (c.f. Eq. [T5l) . 
The external legs are hard and nearly on-shell due to pinch singularities, i.e. k ~ p ~ T and 
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k r a p k r a p 

FIG. 12: Possible assignments of r,a labels for rung (vl); other assignments are subleading because 
of, e.g. the absence of pinching propagators inside the soft vertex correction. Expanding the rr 
propagators as G r ^ F (k) = (1 ± 2n B / F (k )) (G r ^ F (k) — p (k)j , we can show that these two 
contributions nearly cancel (c.f. Eq. (|34h ). Figure taken from [271 ] . 



k r a P 




k r a r a p 

FIG. 13: Momentum and r,a labels used to do the power counting of rungs (u2) (c.f. Fig. [9] and 
Eq. ([35]) ). A thick line corresponds to a soft propagator. 

k 2 ~ p 2 ~ 0(e 2 T 2 ). We consider the case where the loop momentum I is soft (the case where 
I is hard corresponds to the non-collinear rung (e2)). Similarly to rung (ul), the integral 
over dl° is dominated by the region when the two propagators G c ^(k — p + l)G r F (k + I) 
nearly pinch: this happens when the hard photon is collinear with the hard electrons (i.e. 
when ~ 0(g)). To get a leading order contribution, we also need / to be spacelike, 
giving G^ HTL (/) ~ 0(e~ 3 ). Note that G r p{k — p) does not have any pole in 1° and thus 
do not give any singular contribution to the integral. Collecting all powers of e, we get 
e 4 x e 2 x (e 2 x e 3 ) x e -3 x (e -2 x e -2 ) ~ 0(e 4 ), as for rung (ul). Thus rungs (u2)-(v2) must 
be included in the leading order collinear analysis. Note that rungs (w2)-(x2) are already 
included in Sect. IIV Al (compare with rungs (e2) and (f2)). 
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More generally, power counting allows all fermionic self-energies with any number of 
vertex corrections to the hard photon to be included (some examples of self-energies are 
shown in Fig. [T3J along with the corresponding rungs). To see that, note that a vertex 
correction (on any end of the hard photon) adds two explicit vertices (e 2 ), an integral over 
soft momenta (e 3 ), two collinearly singular propagators (e -2 ) and one soft rr propagator 
(e~ 3 ). Thus adding a vertex correction amounts to an 0(e°) change to the self-energy. This 
is true for any number of vertex corrections on both ends of the hard photon. 

However, most of the rungs in Fig. [TD do not contribute at leading order and can be left 
out of the analysis. For instance, using arguments similar to the ones used for rungs (vl)- 
(wl), it is easy to show that cancellations between different r,a assigments occur for rungs 
(ttl), (uul), (yyl), (zzl), (w2) and make them subleading. Rungs (vvl), (wwl), (xxl) are 
also subleading since there is no hard exchange between the side rails (same argument as 
for rungs (yl)-(wl)). Lastly, rungs (tt2), (uu2), (ww2) should not be considered in the first 
place, since they are obtained by opening fermionic self-energies with respect to soft HTL 
resummed spectral densities; this not allowed, since the functional derivative in Eq. (|27j) is 
with respect to ps, not Pbktl- Note that the above arguments are robust with respect to 
the addition of more vertex corrections to the self-energies. Thus after a careful analysis, 
we see that all the rungs in Fig. [TJ]are subleading, except rungs similar to rungs (ul), (xl), 
(u2) and (v2) with an arbitrary number of crossed soft photons. This completes the power 
counting for the class of rungs obtained from fermionic self-energies. 

The second class of rungs comes from opening bosonic self-energies (see Fig. [TO}) . We 
first notice that, as in the non-collinear case, rungs (u3)-(x3) are the "mirror images" of 
rungs (u2)-(v2). Consequently, the power counting of rungs (u3)-(x3) is similar to rungs 
(u2)-(v2) and shows that they are all 0(e 4 ). As for fermionic self-energies, an infinite 
number of bosonic self-energies must also be included at leading order. These self-energies 
are shown in Fig. [15] and are analyzed in details in 58] . The key point is that the addition 
of soft exchange photons to the bosonic self-energies in Fig. [15] are in fact corrections to the 
incoming/ougoing hard photon vertex; they are similar to the multiple vertex corrections in 
the fermionic self-energies (c.f. Fig. |L4"|) and can be treated in the same way. 

Lastly, rungs (u4), (v4) (c.f. Fig. fTOj) are subleading and are not considered in the collinear 
analysis. The power counting for these rungs is very similar to the one for (ul), (xl), with 
the difference that the two bosonic exchange propagators are replaced by two fermionic 
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FIG. 14: Possible 3-loop fermionic self-energies that contribute at leading order. A thick line 
corresponds to a soft propagator. The corresponding rungs (obtained by functionally differentiating 
the self-energies with respect to pf or pg) are also shown. Note that not all the possible rungs 
are shown, this is only a representative sample. Power counting arguments tell us that only rungs 
(ssl) and (ss2) must be included in the collinear analysis. 

propagators. From a power counting point of view, it implies that rungs (u4), (v4) do not 
have the Bose-Einstein enhancement associated with a soft bosonic exchange propagator and 
are thus subleading. In any case, rungs (u4), (v4) are not allowed since they are obtained 
by opening bosonic self-energies with respect to ps htl and not p B as prescribed by Eq. 
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FIG. 15: Possible 3-loop bosonic self-energy that contribute at leading order. A thick line corre- 
sponds to a soft propagator. The corresponding rungs (obtained by functionally differentiating the 
self-energy with respect to pf) are also shown. Note that not all the possible rungs are shown, this 
is only a representative sample. Power counting arguments tell us that only the rungs with soft 
corrections to one vertex at a time must be included in the collinear analysis. 

The infinite number of rung types with arbitrary number of soft vertex corrections must 
be resummed using another integral equation. This is the origin of the comment made in 
Sect. Ill Bl about the necessity of using two integral equations embedded in each other to 
obtain a leading order result. The fact that an arbitrary number of soft vertex corrections 
are allowed is a direct consequence of collinear singularities and is the diagrammatic imple- 



mentation of the LPM effect 
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60] . We will see more clearly in Sect. IV Bl how these 
collinear rungs are related to the 1 + iV — > 2 + N collinear scattering processes found in e.g. 
Ref. 161] . 



V. DERIVATION OF THE INTEGRAL EQUATIONS 

A visual summary of the resummation program for the computation of shear viscosity 
at leading order, including all the necessary collinear and non-collinear rungs, is shown in 
Fig. [T6j The next goal is to write down the appropriate integral equations (including all the 
cuts of the topologies shown in Fig. fT6l) and show their equivalence to the results of Arnold, 



Moore and Yaffe [4], 125] obtained using effective kinetic theory. 

The essential steps leading to the linearized Boltzmann equation for shear viscosity in 
scalar theories can be found in 0|- The same method (generalized to fermions and 
gauge theories) is used in 26j to obtain the Boltzmann equation for electrical conductivity 
in hot QED. Since the details of the method for both bosons and fermions are available in 
the litterature, we only outline the important steps in the present section. 

As is done in Sect. IIVI we separate the analysis of collinear and non-collinear rungs in 
Eqs. ( |T5l) -( |T6l) according to /C = {M. + A/)^ 7 , where Af and M. correspond to the collinear 
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and non-collinear rungs, respectively. 



A. Integral equation without collinear singularity (M) 

We first consider the case with only non-collinear rungs for simplicity. We want to 
write down Eqs. (|T5l) - (jT6l) in the "1-2" formalism, including all possible cuts of the effective 
vertices, the side rails and the rungs. A compact way of doing that, pioneered by Jeon 18], is 
to write the integral equations ([T5l) - ([T6l) in matrix form, where each component corresponds 
to a different cut. We can thus write: 

T> F (k) = I F (k) + J M w (k,p)r F (p)V P (p) 

+ JwT 4 M ( J i (k, ^* (p) (36) 



7> B (k) = X B (k)+ I -0^M (3) (k,p)F F (p)V F (p) 



+ S Mw{Kp)F B (jp)T> B {p) (37) 



where we use boldface letters to emphasize the fact that T> and X are 4 component column 
vectors and 7Q = M-usF are 4x4 matrices. See Fig. [17] for the graphical representation 
of this equation. To avoid cluttering the equations with indices, we do not write Lorentz 
indices explicitly (see Eqs. (1221) - (1231) for the correct assigment of Lorentz indices). We again 
emphasize that, here and in the following, the order of the various fermionic components 
is not respected; for example, the effective vertex TD F (jp) has a Dirac structure and should 
be sandwiched between the two fermionic propagators contained in ?~(p), something that 
is not apparent from the present notation. Only explicit calculations show that the Dirac 
structure all works out. 

The matrix equations (l3T)|) - fl3"TI) are amenable to some simplifications. First note that the 
side rail matrix T can be decomposed into outer products in the following way: 

Fb/f(p) = w B/F (p)u T B/F (p) +h B/F (p)j B/F (p) (38) 

It can be shown that h B / F {p) and j b /f(p) do not contribute to the final integral equation 
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FIG. 16: Diagrammatic summary of our leading order calculation of shear viscosity in hot QED. 
The basic integral equations (c.f. Eqs. f)15[) - (|16[) ) are due to the usual pinch singularities; their 
solutions (represented by a grey half circles) must be substituted in the initial Kubo relation to get 
the shear viscosity. All the rungs included in the kernels of the basic integral equations (represented 
by grey rectangles) are shown. These rungs can be divided in two categories: those corresponding 
to 2 — ► 2 scatterings (located outside the dotted rectangles) and those corresponding to 1 — > 2 
collinear scatterings (located inside the dotted rectangles) . The collinear rungs represent an infinite 
number of vertex corrections (represented by a grey triangle) that are resummed using the collinear 
integral equation (c.f. Eq. (06])). Also shown are the self-energies that must be resummed in the 



side rail propagators in order to satisfy the constraints due to the Ward-like identities (|26 p - ([29 
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FIG. 17: Graphical representation of the integral equations (I36p - (|37p , Only the term /Ci is repre- 
sented here; the three other terms have the same cut structure. The inhomogeneous term T(p) is 
not shown explicitly. It is a column vector with only the first component being nonzero, since it is 
an operator insertion and is thus point-like (i.e. it can only be in the shade or out of the shade). 



Figure taken from Ref. 



26]. 



13,|2a]. The expressions for ub/f(p) an d wb/f(p) are: 



u>b/f(p) 



2n B/F (p°)A+(p) 



p B/F 



1 1 



U 



B/F 



(p) = ±(l 1 ^) 



(39) 



Multiplying Eq. fl36|) by u F (k) and Eq. fl37|) by u^(k) from the LHS, it is possible to re- 
duce the matrix integral equations into one-dimensional integral equations (in "1-2" space). 
After this operation, the reduced rung kernels Ug, F (k)A4^(k,p)wB/F(p) becomes a sum 
of different cuts. Using unitarity and KMS relations for 4-point functions in addition to 



some special relations between the rungs 
simplified. After some algebra, we get: 

d 4 p 



m, 



V F (k) = -I F (k) 



V B {k) 



+ 

-T B {k) + 
+ 



(2tt)< 
d 4 p 

d 4 p 



1 + e" 



261 ]. the reduced rung kernel can be further 



M$(k,p) + e^M£(k,p)) A f (p) 



/3fc° k AX), 



:i + e 



-/3k -' 



M$(k,p) + e^M%(k,p)) A+(p) 



(2tt)< 



(1 - e-e*) (M%(k,p) + e^'Mf>(k,p) ) A+(p) 



,(3) 



2E p 'y p 
V F (p) 



(40) 



d 4 p 
J2n) 4 



:i-e 



-0k^ 



M$(k,p) + e^M%>(k,p)) A+(p) 



f>k° A ,(4), 



2E P 7p 



(41) 



where £>_b/f = Tu%/ f 'Db/f, the JAi/s = 0, ...,3) refer to the matrix components of 
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A4 (correspond to the different ways of cutting the rung kernel, see Fig. IT7I) . we defined 
Ap(p) = fiAp(p) and we used the fact that Up oc 7*. 

Equations fHUl) - f|41l) . with the non-collinear rungs in Fig. [TH] and the current insertions 
given by Eq. ([TBI) , contain the necessary ingredients for computing the shear viscosity at 
leading order (neglecting collinear physics). One could in principle solve these equations 
numerically to obtain the effective vertices and then compute the shear viscosity from the 
Kubo relation ffT2]) . 

In this paper, we instead want to show the equivalence between Eqs. (T40l) - fHTl) and the 
effective kinetic equations of Arnold, Moore and Yaffe 4, 2s|. To do that, we first note 
that the non-collinear rungs in Fig. [TH] are all made of 4-point functions with two external 
vertices in the shade and two out of the shade. Since cut propagators represent nearly on- 
shell thermal quasi-particles, the non-collinear rungs can be naturally interpreted as 2 — ► 2 
scattering processes. A straightforward calculation shows that the sum of all the leading 
order non-collinear rungs contained in the kernels -M(i) + M.(2) and A4(3) + M.U) can be 
converted into the square of a scattering matrix, where the scattering processes are given in 
Figs. [T8lfT9l It is easy to see that diagrammatically, starting from the scattering processes 
(see the captions of Figs. [T8lfT9l for details). Multiplying Eq. [40] from the left (right) by the 
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spinor u x (k) (u x (k)), i 
D e (k) = -I e (k) + 

x 
+ 

X 

D 7 (fc) = -I 7 (A;) + 
x 
+ 

X 

where the M uvmn (k,p; I, /')'s are the 2 — > 2 scattering processes shown in Figs. fT8lfT9l we have 
used explicit labels for excitations (e =electron, p =positron, 7 =photon), we have defined 
Dp{k) = u x (k)Vp(k)u x (k) and redefined Ds(k) = T>B{k). The sums are over flavors (f), 
species (s) and helicities (h). 

We also note that the integrals in Eqs. (l40l) - ()4Tl) are over d 4 p, while kinetic equations 
deal with on-shell quasi-particle propagation. Thus to make the connection between the 
two approaches clearer, we use the delta functions present in the cut propagators to put 
excitations on-shell in Eqs. fj40l) - (|4TI) . We remark here that at finite temperature, both 
positive and negative energies appear in the cut lines. In addition, the distribution functions 
inside the cut propagators provide the thermal phase space for the gain and loss terms in 
the kinetic equations. 

Defining the deviations from equilibrium XB^k) = Dsik) flE^k and XF^k) = Dp(k)/rk 
and noting that 7^ and 1^ can be written as squares of scattering processes, we follow the 



e get: 

d 4 p dH dH' 

(2tt) 4 (2vr) 4 (2tt) 4 

f,s,h 



(27r)V 4) (fc + p-/-/')(l+e 



E 

v£(e,p) ; m,ne(e,p,7) 



M evmn (k,p;l,l')\ 2 A+(-p)A+(/)A+(0 



D v (p) 



1 f,s,h 

- \M eimn (k,p-JJ')\ 2 A+(-p)A+(J)A+(J' 



m,n 



2E P 7p 



(42) 



M, vmn (k,p-JJ')\ 2 A+(-p)A+(i)A+(0 



f,s,h 

E 

■uG(e,;p) ; m,ng(e,p,7) 

(&(Si(&w vi, (* + "-'- i ')( i - 

- |M 77mn (fc,p;/,/')| 2 A+(-p)A+(0A+(O 



A,(p) 



gr(g) 

2£"p7p 



(43) 
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FIG. 18: Collision term of the (fermionic) linearized Boltzmann equation including only leading 



order 2 — > 2 scatterings 



251 ] . The labels it-n represent the species (fermion or photon), flavor and 



helicity of the excitations (here u = electron) . A straightforward calculation shows the equivalence 
between the scattering processes and the rungs of Figs. 0H5J The correspondence goes as follows 
(the letters refer to diagrams in Figs. H1I51 and [TBI respectively): (al) = (b) 2 and (d) 2 , (bl) = (e) 2 
and (f) 2 , (cl) = (a) 2 and (c) 2 , (dl) = (g) 2 and (h) 2 ,(gl) = (g)(h) and (h)(g), (jl) = (e)(f) and 
(f)(e), (kl) = (b)(a) and (d)(c), (11) = (a)(b) and (c)(d), (c2) = (j) 2 and (l) 2 , (d2) = (i) 2 and (k) 2 , 
(e2) = (i)(j) and (k)(l), (f2) = (j)(i) and (l)(k). 
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FIG. 19: Collision term of the (bosonic) linearized Boltzmann equation including only leading 



order 2 — > 2 scatterings 4, 25]. The labels u-n represent the species (fermion or photon), flavor and 
helicity of the excitations (here u = photon). A straightforward calculation shows the equivalence 
between the scattering processes and the rungs of Figs. [7)0 The correspondence goes as follows 
(the letters refer to diagrams in Figs. I7H81 and \W\ respectively): (a3) = (e) 2 and (g) 2 , (d3) = (b) 2 
and (d) 2 , (e3) = (a) 2 and (c) 2 , (h3) = (f) 2 and (h) 2 , (13) = (e)(f) and (g)(h), (k3) = (b)(a) and 
(d)(c), (m3) = (a)(b) and (c)(d), (o3) = (f)(e) and (h)(g), (a4) = (i) 2 and (l) 2 , (b4) = (j) 2 and 



(k) 2 , (c4) = (i)(j) and (l)(k), (d4) = (j)(i) and (k)(l). 
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same procedure as in Ref. [26J] to get: 

d 3 p d 3 l d 3 l' 



-(l-n e (k ))I e (k) = - I , d J , d ! , d J (27r)V 4) (fc + p-/-n 
v e\ }) e\ ) 2 j ( 27 r)32£ p (2n) 3 2Ei (2tt) 3 2E v v ; v y ' 

f,s,h 

x |M e , mn (A;,p;/,/ / )| 2 n e (fcV(j9 )(l±n m (/ ))(l±n n (/ /0 )) 

v,m,n 

x [*.(*) + Xv (p) - Xm (l) - xJl')] (44) 
-<1 WOW*) = \ J ( $ 2Er {2 ^ 2El ( J4 f + "-'-'') 

x ^ |M 7 „ m 4fc,p;/,r)| 2 n 7 (fc )n,(p )(l±n m (/°))(l±n n (/' )) 

v,m,n 

X [X 7 ( fc ) + X«(P) - Xm(l) ~ Xn(l')] (45) 

These coupled integral equations are identical to the ones obtained by Arnold, Moore and 



Yaffe 



251 ] using effective kinetic theory (without collinear processes). 



B. Integral equation with collinear singularities (AT) 

The case with collinear rungs can be treated using the same tools used for non-collinear 
rungs. Since the steps are essentially similar, we focus in this section on how to write the 
collinear version of Eqs. (H2]) - ff^3|) and how to properly write the leading order collinear rungs 
asl + iV^2 + iV r scattering processes. 

According to the power counting of Sect. IIVBI the rung kernels are given by A/(i) = 
(til) + (xl), jV(2) = (u2) + (v2) and jV(3) = (u3) + (t)3) + (w3) + (x3), along with all similar 
rungs with an arbitrary number of soft photon corrections (see Fig. |2"U|) [64| . For Nm and 
A/(3), this infinite number of soft photons corresponds to vertex corrections. This is also 
the case for J\f(x), although it is less apparent; the key point to understand this is to note 
that the reduction of the side rail matrix JF into a cut propagator in the pinch limit (c.f. 
Eqs. ( l38l) - (!39l) ) effectively closes the rung, making the soft exchange photon between the two 
side rails look like a vertex correction (see Fig. I20i) . 

This infinite number of vertex corrections due to collinear singularities can be resummed 
by defining an effective vertex VV(/c,p) that satisfies the following integral equation: 

V F {k,p) = I F {k,p) + J j^^M C o\\{k,p,q)J r F{k,p 1 q)V F {p,q) (46) 
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FIG. 20: Rungs that are part of the leading order collinear analysis. A thick line corresponds to 
a soft propagator. Due to collinear singularities, adding a soft photon is an O(l) correction; thus 
for each rung kernel , an infinite number of rung types contribute at leading order. These soft 
vertex corrections can be resummed using another integral equation (c.f. Eq. (|46|) ). For J\fm, it is 
not readily apparent that these soft photons correspond to vertex corrections. To see that, notice 
that the rung is effectively closed with a cut propagator coming from the reduction of the side rail 
matrix T (c.f. Eqs. ([55]) - using this fact, we see that the resummed JVm rung is equivalent 
to the one in parenthesis. 



where A/" co ii is a rung with a single soft photon exchange and the external photon is collinear 
with the electron. See Fig. [2T] for an illustration of this integral equation. This resummation 
is relevant for photon production including the LPM effect and is done in great detail in 



58, 



621 ] . Instead of using the closed-time-path or "1-2" formalism, they use the r/a formalism 



and are able to put the integral equation in a form convenient for numerical purposes. 
Applying the reduction procedure of Sect. IV Al to the collinear rung kernels of Fig. [20, 
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FIG. 21: Illustration of the integral equation (|46p . The kernel Af co n is made of only a single soft 
photon exchange. 

it is possible to write a set of integral equations similar to Eqs. ( )40l) - (l411 (with the A4's 
replaced by A/"'s). Note that A/32 = A/32 = (that particular cut does not exist for 3-point 
functions) and A/32 is subleading; there is thus no A/32 contribution in the collinear version 
of Eqs. gOD-flUD. 

To obtain the collinear version of Eqs. P2|) - fj43|) . we need to express the collinear rungs 
in Fig. [20] in terms of 1 — > 2 collinear scattering processes. We first note that all the 
collinear rung kernels can be written as the multiplication of a bare 1 — > 2 vertex V^, 2 
and a resummed vertex V F . For instance, we have A/^( ) = i~ ^1-^2] nvft ~~ P) [~ ^f\* 
and A£j2 (5) = [ — ^f\ ~~ P) l~~ ^1-^2]* ■ Note also that A/^ ^ and A/^* ™ are complex 

conjugate of each other; this is the consequence of the trace over the fermion structure in 
the Kubo formula, which allows us in this case to freely interchange the Dirac structure. 
The sum of A^L) an d ■N'^tb) * s ^ nus twice the real part of N^, a y The other rung kernels 
AA 2 ) and J\f^ can be similarly written. Following the procedure in Sect. IV A\ we obtain: 

d 4 p dH 



D e (k) = -I e (k) + 



(2vr) 4 (2tt) 4 
x [2Re[V£V^ 2fi r n (k;p,l) A+(p)A+(Z 
d 4 P dH 



(2tt) 4 (2%y 



(27r) 4 5 (4) (A;-p-/)(l+e 



D^k) = -I 7 (k) + 



[2He[Vg^ 2 J%(*;p l A+(p)A+(J)] 
(2vr)V 4 )(A;-p-/)(l -e 



d 4 p dH 



(2vr) 4 (2tt) 4 



ge(g) 

2E p j p 



(47) 



£ 2Re[V^ 2 J^(^;P,0 A+(p)A+(Z) 



AXp) 



(48) 



Figure 
Fig. 



shows the 1 — > 2 collinear scattering processes corresponding to the rungs in 
These scattering processes are the ones that are responsible for the LPM effect 
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6l| . The equivalence with kinetic theory can again be seen by defining the deviations from 
equilibrium Xb(^) = -Ds(^)/2-E P 7fc and Xf(^) = Dp(k)/Tk and following the procedure in 
Sect. IV Al The final result is: 

-(l-n e (k°))I e (k) = \J ^ { Jy 2Ei P*)W\k - p - I) £ 2Re[V^ 2 JU*!?. 

x [n e (k°)(l ± n v (p°))(l ± n m (/ )) ( Xe (k) - Xv (p) ~ Xm(0)] 

x [n e (k°)n v (p°)(l ± n m (/ )) ( Xe (fc) + x„(p) - Xm(0)] ( 4 9) 
-(l + n 7 (0)W = \J { J^ P 2Ep { J** 2Ei (2n)W\k -p-l) f2^e{V^ 2 ,)Uk-,P,l) 



x [ nj (k )(l ± n v (p°))(l ± n m (Z )) ( X #) - Xv (p) - Xm (Z))] 



+ 



rf 3 p dH ^ 



(2^(2^ (^(* + ^-0 E^e[V^ 2 ^(^;0 
x [n 7 (A;V(p )(i ± n m (/ )) + X „(p) - Xm (Z))] (50) 

This last equation is identical to the linearized Boltzmann equation with the transverse 



momenta non-integrated of Arnold, Moore and Yaffe 



25] obtained using kinetic theory. 



VI. CONCLUSION 



In this paper, we have derived the integral equations needed for the calculation of shear 
viscosity at leading order in hot QED. A visual summary of our calculation is presented in 
Fig. [THl Our calculation includes all the necessary 2^2 scatterings and 1 — >• 2 collinear 
scatterings responsible for the LPM effect. The equivalence with the effective kinetic theory 
results of Arnold, Moore and Yaffe has also been established. 

We used Ward-like identities to put constraints on the form of the rung kernels. In 
the case of electrical conductivity, these constraints are directly related to gauge symmetry 
and gauge invariance is explicitly enforced. The issue is not so clear in the case of shear 
viscosity, since the Ward-like identities used to obtain the constraints derive from space-time 
symmetries. But since the constraints (I26I) - (I29I) have the same form as those obtained from 
the usual Ward identity in the case of electrical conductivity, then Eqs. (126]) - (129]) should be 
sufficient to preserve gauge invariance; we have shown explicitly that this is indeed the case. 
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c.c 
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FIG. 22: 1 + N — > 2 + N collinear scatterings included in Eqs. P7]) - P5]) . Scatterings with soft 
photons are represented by a wiggly line ending with a crossed circle. These collinear processes 
correspond to the ones identified in Ref. [6l|] as responsible for the LPM effect. 

The formalism developped in the present paper could in principle be applied to other 
transport coefficients, such as the bulk viscosity. In scalar theories, bulk viscosity is sensi- 
tive to number changing processes at leading order whereas shear viscosity is not, making the 
calculation of the former more involved. In gauge theories, all transport coefficients are sen- 
sitive to these processes even at leading order. Thus only minor modifications (e.g. changing 
the current insertion) are required to adapt the present formalism to the computation of 
other transport coefficients in QED. 

Other extensions include the generalization of our formalism to non-abelian gauge theories 
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and the computation of transport coefficients beyond leading order. Of course, generalizing 
our calculation to say QCD would be quite demanding, because of the presence of ghosts 
and a more complicated collinear resummation (due to the self-interaction of gluons among 
themselves). Going beyond leading order would also be very challenging, if just for the very 
subtle power counting that would be required or the possibility of "new physics" (e.g. new 
singularities that would require additional resummations) appearing in the calculations. 
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n 

Note that rungs (vl), (wl) were included in [26J instead of rungs (ul), (xl). As explained 
in Sect. IIVB) the size of both sets of rungs is (naively) 0(e 4 ), but (vl),(wl) turn out to be 
subleading due to a cancellation. However, since both sets of rungs can be seen as vertex 
corrections, the final result is the same f^ . 
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APPENDIX A: GAUGE INVARIANCE OF SHEAR VISCOSITY 



In this appendix, we present arguments showing the independence of the shear viscosity on 
the gauge parameter £ at leading order in QED. The starting point is the Kubo relation (|14p 
with the accompanying integral equations (Tl5l) - ([T6"]) . Since the current insertions Ib/f do 
not depend on £, the goal is to show that the side rails Tbif arid the effective vertices 
V B /p do not have any gauge parameter dependence. The side rails Tp may depend on £ 
through their resummed self-energies (we come back to Tb later). The effective vertices are 
solution to the integral equations f|T5l) - f|T6l) ; V b /f may thus depend on £ through the side 
rails J-'b/f and the four rung kernels M.. Since the Ward-like identity constraints fT26T) - (l2§|) 
establish a direct relation between the rungs and the self-energies resummed into the side 
rail propagators, the only thing we have to show is that the resummed self-energies are gauge 
parameter independent. 

Before presenting our arguments, note that there are two different types of Ward identities 
in this study. The first one is k^A^ = and is discussed in e.g. 511] ■ The second one is 
kf/D^ oc T,j(k) and is discussed in Sect. IIHI The latter is valid only in the kinematic limit 
relevant for transport coefficients and when the two external resummed propagators pinch. 

For definiteness, let's consider fermionic self-energies (bosonic self-energies can be ana- 
lyzed in a similar way). The leading order fermionic self-energies included in the side rail 
propagators are shown in Fig. [16] They can be divided into two subgroups, those containing 
and those not containing collinear singularities. Using the usual Ward identity k^A^ = 
(where A 11 is any amplitude with all external excitations on-shell), we can show that all the 
non-collinear self-energies are gauge parameter independent. For example, the first fermionic 
self-energy in Fig. [16] contains two uncut (i.e. off-shell) photons. The £ dependent part in 
the photon propagator is proportional to t^k^k 1 * and we can send each of those fe's to a ver- 
tex. Since all the electrons are cut (i.e. on-shell), the four vertices dotted with a k are zero 
due to the Ward identity. The first self-energy is thus £ independent. The second fermionic 
self-energy in Fig. [16] contains two cut photons that correspond to external (thermal) exci- 
tations. From the Ward identity, any amplitude dotted with a k^ is zero; thus the ^k^k u 
part of the cut propagators does not contribute. Using similar arguments, we can show that 
the other two non-collinear fermionic self-energies are also gauge parameter independent. 

The case of collinear fermionic self-energies is slightly more involved than the non-collinear 
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one, since self-energies with an infinite number of soft corrections to the hard photon vertex 
contribute at leading order. To deal with this infinite number of soft corrections, note that 
all the electron propagators are nearly on-shell (i.e. on-shell to within 0(e 2 T 2 )), since the 
hard photon is collinear to the incoming electron and the soft photons can only push it 
out of his mass shell with 0(eT) kicks. This implies that the electron wavefunction obeys 
a "nearly on-shell Dirac equation" (J> — m)u(p) ~ 0(eT). From this nearly on-shell Dirac 
equation, we deduce that the Ward identity obeyed by each vertex with two incoming nearly 
on-shell electrons is k^A^ ~ 0(eT). As for the non-collinear case, we can use this "nearly 
on-shell" Ward identity to deal with the ^k^k 1 ' term in each photon. Since we must use 
the nearly on-shell Ward identity twice for each photon, the £ dependent part of the self- 
energy is suppressed by at least 0(g 2 ) compared to the rest of the diagram. Thus the gauge 
non-invariance of collinear fermionic self-energies only appears at higher orders. 

Another potential source of £ dependence is Tb (he. the bosonic side rail propagators). 
To see that Tb does not give rise to any £ dependence, note that two pinching propagators 
can always be expressed as a cut propagator (c.f. Eq. (1391)). Consequently, the ^k^k u term 
in Tb is always on the external leg of an amplitude and vanish by the Ward identity. This 
completes our argumentation and shows that the shear viscosity is gauge invariant at leading 
order. 
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